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Abstract -We extend the theory of quasipotentials in dynamical systems by calculating, within 
a broad class of period-doubling maps, an exact potential for the critical fluctuations of pitchfork 
bifurcations in the weak noise limit. These far-from-equilibrium fluctuations are described by finite- 
size mean field theory, placing their static properties in the same universality class as the Ising 
model on a complete graph. We demonstrate that the effective system size of noisy period-doubling 
bifurcations exhibits universal scaling behavior along period-doubling routes to chaos. 


Period-doubling bifurcations have been observed in a wide variety of natural systems 
spanning many areas of science 1 . Univariate, discrete-time maps are the simplest dynam¬ 
ical systems to exhibit a period-doubling route to chaos [2j[3j. Applications range from the 
population dynamics of species with non-overlapping generations [IJ[5] to the oscillations of 
rf-driven Josephson junctions [6,7|. The impact of noise on period-doubling maps has been 
extensively studied in both ecology [8 11 and physics 12 30 , including the universal scal¬ 
ing of the Lyapunov exponent along period-doubling routes to chaos 31-36 . Connections 


between such noisy dynamical systems and the universality classes of equilibrium statistical 
physics has been a subject of great fascination 37 -49 . The theory of quasipotentials, pro¬ 


viding a formal link between equilibrium and nonequilibrium physics, has been applied to 
In particular, the theory of quasipotentials has been used to estimate 

and to estimate the 


many systems 


50 51 


escape times from the attractors of noisy period-doubling maps 52-54 
invariant probability distributions of their strange attractors in the chaotic regime 53|55|56 


Here, we investigate the invariant probability distributions that characterize critical fluc¬ 
tuations in the pitchfork bifurcations of period-doubling maps far from the chaotic threshold. 
In the limit of weak noise, we find an exact correspondence between the static behavior of 
fluctuations at a pitchfork bifurcation and the critical behavior of finite-size mean field the¬ 
ory 57 . This correspondence places pitchfork bifurcations in the same universality class 
as the Ising model on a complete graph 58 -62 . Analytical estimates of critical exponents 


and amplitudes agree well with the results of numerical simulations. We conclude with evi¬ 
dence of universal scaling behavior in the effective system size of critical fluctuations along 
period-doubling routes to chaos. 
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We start by considering a noisy, one-dimensional map that can be written as 


x t+ i = + Xa(r,x t )£,t, 


(1) 


where the a; t are random variables, r is the control parameter, /.t(r, Xt ) is a deterministic 
period-doubling map, A > 0 is an overall constant factor setting the noise level, a(r,x t ) is a 
nonnegative function, and the are independent, identically distributed standard normal 
random variables. We restrict our attention to fluctuations arising from the effects of weak 
noise on supercritical pitchfork bifurcations. A pitchfork bifurcation in a period-doubling 
sequence corresponds to a loss in stability of the steady state fixed points, x*, with ( r c ,x *) 
being the critical point where [l 63,64] 


M(r c ,x*) = x*, 

d x /x(r c ,x*) = -1. 


( 2 ) 


Note that x* is the stable fixed point for r < r c , and we define it to be the unstable fixed 
point for r > r c . The first bifurcation breaks a time-translation symmetry: below r c , the 
stable attractor is a steady-state fixed point, and trajectories on the attractor are invariant 
to all time-translations; above r c , the stable attractor is a two-cycle and corresponding 
trajectories are only invariant to translations by an even-number of time steps. We will 
investigate the critical fluctuations that emerge from this Zi symmetry breaking in the 
limit of weak noise (A —> 0). Our results rely on the additional assumptions that [i(r,Xt) is 
a C 3 


map 


65 


i.e. three-times differentiable with respect to Xt, and that the stable basin of 
attraction of /z(r, x t ) is bounded on some finite interval. 

Noisy maps in the form of Eq. |T] have been applied to the study of stochastic fluctua¬ 
tions in ecological population densities 8,10 . If p(r, Xt)=rxt(l — Xt) and a(r,Xt)= 1, then 
Eq. [I] is the logistic map with extrinsic additive noise. Non-constant choices for a(r,x t ) 
model intrinsic noise or a combination of intrinsic and extrinsic noise. When applied to 
population dynamics, common choices might be cr(r,Xt)=n(r,Xt), modeling multiplicative 
environmental noise, or er(r, Xt.)=y/ Xt), modeling demographic noise. 

Over long enough time scales, a noisy trajectory, x t , will depart from the stable basin of 
attraction of n(r,Xt). However, for small enough A and for initial conditions on the interval 
bounding the stable attractor of n(r,Xt), a metastable probability distribution for Xt will 
persist over “intermediate” time periods that are much longer than the duration of initial 
transients but exponentially shorter than the characteristic time to depart from the stable 
attractor. 

This separation of time scales allows us to quantify the critical fluctuations of pitchfork 
bifurcations. The metastable distribution of the map variable, Xt, is approximately sta¬ 
tionary at any given intermediate time, T. We define the pitchfork order parameter, m, 
as 

(3) 


m = ^Ot+i - x T ), 


where a is a model-dependent scale factor that maps the pitchfork order parameter onto the 
Ising mean-field magnetization. To obtain numerical estimates of the order parameter, we 
use ergodicity 55 and sample statistics along a single trajectory of Eq. [l] at intermediate 
time scales. 

Fig. a is a histogram of the unsealed order parameter, m/a, for a logistic map with 
weak additive noise. The distribution of the order parameter is highly symmetric and can 
be well-approximated by a pure quartic. These numerical results motivate us to seek an 
analytical correspondence between the critical fluctuations of pitchfork bifurcations and 
those of finite-size mean field theory for the Ising model. The latter are described by the 
distribution 

Pmft(to) oc exp , (4) 
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Fig. 1: (Color online) Histogram of the unsealed order parameter, m/a, at a noisy supercritical 
pitchfork bifurcation. Blue points and error bars are estimated from numerical simulations of Eq. [l] 
at the first period-doubling transition of a logistic map, p(r, x)=rx(l—x) and r=r c = 3, with additive 
noise, A=10~ 4 and a{r, x)=l. Similar results hold for all A < 10 -2 . The black line is a fit of a 
pure quartic distribution to the simulation results (chi-square per d.o.f. 1.3). Fluctuations in the 
pitchfork order parameter described by a quartic distribution suggest a correspondence between the 
critical behavior of noisy pitchfork bifurcations and finite-size mean field theory. 


where N is the number of spins and the free energy (the large-deviation potential) at reduced 
temperature t is 

f(m) = ^tm 2 + ^m 4 , (5) 

to leading order in N. 

The first step is to calculate the metastable distribution, h(x), of the map variable, x, 
near the first period-doubling pitchfork bifurcation in the weak noise limit. We re-write the 
noisy map of Eq. [I] in integral form and obtain a self-consistent equation for h(x) 


h(x) = J dx'h(x')k(x',x), 

where the kernel describing the behavior of the noisy map is given by 


k(x',x) = 


\/27rA er(r, x') 


exp 


(x - M t , x')f 
2A 2 cr(r, x') 2 


( 6 ) 

(7) 


For weak noise, h{x) will be sharply peaked near x*, and we can take the limits of integration 
in Eq. [6] from — oo to oo. Near the critical point, x* is approximately the same as x* with 
the leading order correction given by 


x — x^ = — 


d r li(r c ,x *) 

2 -{- d r d x [ a( vci x c ) 


(r c - r). 


( 8 ) 


We define h(x — x*) = h{x) and make an exponential ansatz for h given by 

h{z) oc exp (- 1 - [b 2 ,iz 2 6 + b 3A z 3 S + b^ 0 z A + b 5fi z 5 -] - [a\^z H-] ) 

V za a(r c ,x c ) j 


(9) 
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where S = r c — r and the dots indicate higher order terms. In the limit of small A, we apply 
Laplace’s method, and after some algebra, find a self-consistent solution to Eq. [6] where 


^1,0 — *^c) ’ 

62.1 = -Tp(r c ,x*), 

63.1 = £ T ^r c ,x* c )d xx ii{r c ,x* c ), 

Sfi(r c ,x *), 

&5,o = -zSp{r c ,x* c )d xx p(r c ,x* c ). 

6 


( 10 ) 


The expressions in Eq. 10 depend on the critical values of the Schwarzian derivative 65 


Sp,(r c ,x* c ) = - ( d xxx p{r c ,x* c ) + ^ {d xx p(r c ,x* c )) 2 


and the “T-derivative” 


Tn{r c ,x *) = - d r d x p(r, p(r, x))\ ,. 


( 11 ) 


( 12 ) 


The Schwarzian derivative arises from a second-order stability analysis of the normal form 
of a pitchfork bifurcation and is always negative for the supercritical pitchfork bifurcations 
of C 3 maps 65 . The T-derivative arises in our analysis due to the expansion in 5 and is 


always negative for pitchfork bifurcations with a symmetric phase below r c and a broken 


symmetry above, as can be seen from a calculation using the normal form 63 . 

Note that the self-consistent solution for h, as given in Eq. [ 9 ], is not symmetric under 
z —> —z due to the odd terms in z, whereas the order parameter is expected to be symmetric 
under m —> —m, at least to leading order. Given the definition of the order parameter in 
Eq. [3j we have 


p(m) = J J dx! dx 5 [m — — [x — h(x')k(x', x). 


(13) 


Upon calculating the integral, we find that p(m) = pmft(w), as defined in Eqs. [4] and [5] to 
leading order in m and N, with order parameter rescaling factor 


l2Sp(r 


Tp(r c ,x*) 


effective reduced temperature 
and effective system size 

N = 


t = r c — r, 


2A 2 cr 2 (r c , x*) 


(' Tp{r c ,x* c )f 
- Sn(r c ,x*) 


(14) 


(15) 


(16) 


Consistent with previous calculations of quasipotentials in noisy dynamical systems 50 53 
56 , the expression for the effective system size parameter, N in the exponent of Eq. |4| is 


inversely proportional to the variance of the noise, A 2 er 2 (r, x*). Note that the static behavior 
of critical fluctuations at a pitchfork bifurcation is entirely determined by N. In applications 
to ecological systems with demographic noise due to the well-mixed interactions of a finite, 
rather than an infinite, number of individuals 10 , N is directly proportional to population 
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size. Our results quantify the scale of fluctuations for finite-size populations undergoing a 
period-doubling transition. 

For pitchfork bifurcations, the leading odd correction to the potential in Eq.[5|is generally 
0(m7). However, an important quadratic map in ecology is the Ricker map 66 , where 

V{r,x t ) = x t exp(r(l — x t )), (17) 

and this map is symmetric about the first period-doubling bifurcation at ( r c ,x *) = (2,1), 
such that the second derivative d xx n(r c , x*) vanishes and the odd correction to the potential 
at 0(m 7 ) is zero. We speculate that all odd corrections vanish in the potential describing 
critical fluctuations near the first period-doubling bifurcation of the Ricker map. 



jV 1 / 4 


Fig. 2: Finite-size scaling of (a) the absolute value of the pitchfork order parameter, (|m|), and 
(b) the susceptibility, y. Numerical estimates (symbols and error bars) are based on samples 
of the order parameter along Monte Carlo trajectories for the the first period-doubling pitchfork 
bifurcation of a logistic map with weak additive noise, i.e. Eq. [ljwith /i(r,x)=rx(l — x), a(r, x)=l, 
and r=r c =3. Numerical estimates of (|m|) and y are consistent (chi-square per d.o.f. of 0.35 and 2.5, 
respectively) with the predictions of finite-size mean held theory (Eqs. |A2a and A2b| respectively) 
and N as defined by Eq. [16] 
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Fig. 3: Scaling collapse of (a) the absolute value of the pitchfork order parameter, (|m|) and (b) the 
susceptibility, y. Numerical estimates (symbols and error bars) are based on samples of the order 
parameter along Monte Carlo trajectories for the the first period-doubling pitchfork bifurcation of 
a logistic map with weak additive noise, i.e. Eq. [I] with p(r, x)=rx(l — x ), a(r,x)=l, and r=r c = 3. 
Numerical estimates of (|m|) and x are consistent with the predictions of finite-size mean field theory 
(black lines). In panel (b), the vertical grey line marks the peak value of scaled susceptibility, which 
defines the finite-size critical point, r c (N) (Eq. 181. 


We compare mean-field predictions of critical exponents and amplitudes to Monte Carlo 
simulations 67-69 of the first period-doubling bifurcation in a logistic map with additive 
noise. Finite-size scaling (Fig. [2]) and scaling collapse (Fig. [3]) are as predicted by finite-size 
mean field theory. In the finite-size scaling of the order parameter and the susceptibility, as 
defined in Eq. |A1[ numerical results are consistent with mean-field predictions for critical 
exponents and amplitudes. In the scaling collapse of susceptibility, the peak value defines 
the finite-size critical point 

. , 0.88033... 

r c (N) = r c + ■ (18) 


The scaling collapse of Monte Carlo estimates of x* — x* is consistent with the prediction 
of Eq. [8] near r=r c (Fig. [4]). 
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( N 1 ' 4 ) 2 x Reduced control parameter, t 

Fig. 4: Scaling collapse of finite-size corrections to x*. Numerical estimates (symbols and error bars) 
are based on sampling values of x* along Monte Carlo trajectories for the the first period-doubling 
pitchfork bifurcation of a logistic map with weak additive noise, i.e. Eq. [I] with fj,(r, x)=rx(l — x), 
a(r,x)=l, r=r c = 3, and x*=2/3. Those numerical estimates are consistent with our theoretical 
estimate (the first moment of the distribution in Eq. [ 9 ] black line) near the critical point, r = r c - 


We now extend the above results to show that finite-size mean field theory describes 
the critical fluctuations of all the pitchfork bifurcations of a period-doubling map in the 
limit of weak noise. Let r n (dropping the subscript “c” for convenience) denote the nth 
critical point in a period-doubling sequence proceeding from a g-cycle of a deterministic 
map. On r„_i<r<r„, there exist 2 n ~ 1 q stable fixed points of the n-times iterated map, 
Xt+i=H n (r, Xt). At r=r n , stability is lost with pitchfork bifurcations occurring at the points 
x* p , where p&{ 1,..., 2 n ~ 1 q}. With the addition of weak noise, the n-times iterated map 
is, to leading order in A 

x t+ i = Hn(r,x t ) + \<Jn(r,x t )£t, (19) 

where a n (r,Xt) is given by a recursion relation [34] 

cr 2 +1 (r,x) = a n (r,x) 2 (d x p, n (r,fj, n (r,x))) 2 + a n (r,/j, n (r,x)) 2 . ( 20 ) 


Upon substituting Eq. 19 for Eq. [I] and expanding around (r n , x* p ) rather than (r c , x*). 
the above results for the critical fluctuations of the first pitchfork bifurcation in a period¬ 
doubling map apply to all the pitchfork bifurcations in that map. Note that the conditions of 
Eq. [2|apply to all supercritical bifurcations, and that if S/i(r c ,x*)< 0, then Sfi n (r n ,x^ p )<0 


for all n and p 65 


We conclude with a numerical investigation of universal scaling behavior as the first 
sequence of period-doubling bifurcations (g=l) approaches the chaotic threshold at r=roo. 
First, in two different parameterizations of the logistic map, we find evidence (Fig. [5^,) that, 
for large n 


Ifln (j'r 


0 ) oc (r-oo - r n ) oc 8 


Fi 


( 21 ) 


where the Feigenbaum constant <5p=4.6692 ... and the value of —Tp n (r n , 2 * ) is indepen¬ 
dent of p. Second, we investigate scaling in the effective system size of period-doubling 
bifurcations 


N ' l ’ P 2A 2 cr2(r n ,x* ) 


(- Tp n {r n ,x * np )) 2 
-Sfj, n (r n ,x*) 


( 22 ) 
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Fig. 5: (Color online) Numerical evidence of universality in the scaling behavior of (a) 
—T/r n (r n , x„ tP ), and (b) 2A 2 N n . Panel (a) plots values of —T^ n (r n ,Xn :P ) along the first period¬ 
doubling sequence (9=1) in two different parameterizations of the logistic map. Panel (b) plots 
values of 2A 2 N n based on two types of deterministic map, each with three different types of noise. 
In both panels, we find evidence of universal scaling behavior in the the approach to the accumu¬ 
lation point, r=r ao . 


For two different deterministic maps and three different types of noise, we calculate the 
geometric mean of N n p at fixed n 


1/2" 


N n = ( I] N n,P 

,P= 1 


and find evidence (Fig. [5}o) that, roughly 


N n oc 


1 


2A 2 3™’ 


(23) 


(24) 


independent of the details of both the type of deterministic map and the type of noise. This 
result is consistent with the intuition that the noise level, A, must decrease exponentially 
with n to maintain the same effective system size, iV„ iP , at each pitchfork bifurcation. 

We have shown that the critical fluctuations of all noisy pitchfork bifurcations are de¬ 
scribed by finite-size mean field theory for an equilibrium system with N Ising spins. Our 
results were obtained for the case of weak Gaussian white noise but are expected to hold 
more generally for weak noise with a finite second moment and finite correlation time. This 
correspondence confirms the breaking of a Z 2 symmetry at noisy pitchfork bifurcations and 
introduces N as an effective system size parameter that exhibits universal scaling behavior 
near the onset of chaos. In applications to population dynamics with demographic noise, 
N is directly proportional to the population size. Further work is needed to obtain a more 
precise understanding of universal scaling behavior in the critical fluctuations along period¬ 
doubling routes to chaos. 


* * * 
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Appendix A: Summary of results from finite-size Ising mean field theory. 

The existence of finite-size scaling behavior follows from the large-deviation function, Eq. [5] 
of the order parameter distribution, Eq. [4j We will compare numerical estimates for the 
finite-size scaling behavior of noisy pitchfork bifurcations, based on standard Monte Carlo 
techniques 67 69 , to mean-field predictions for the absolute value of the order parameter, 
(|m|), and the susceptibility, %, given by 


x = N ((M 2 ) - < H } 2 ) 


(Al) 


When evaluated at the critical point, t = 0, the mean-field observables depend only on the 
system size, N, namely 

(H)| c 

xlc 

Near the critical point, expressions for observables can be calculated in closed form from 
the moments of Eq. [4] and obey the finite-size scaling relationships 

(M) = (iV 1 / 4 ) _/3 / 1/ G m ^(lV 1 / 4 )i/ I '/j , 

X = (A 1/4 )^G x ((A 1 / 4 ) 1 /"?) , (A3) 

where v = 1/2, /3 = 1/2, and 7 = 1 are the mean-field critical exponents, and the Gx are 
universal scaling functions. We find that Eqs. |A2| and |A3[ with N as defined by Eq. [T6] or 
[22] describe the critical fluctuations of pitchfork bifurcations of period-doubling maps in the 
weak noise limit. 


1 /24tt 5 


r(!) V N 
rm 2 V V2J 


(A2a) 

(A2b) 
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